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Abstract. We analyze the dynamics of a simple but nontrivial 
classical Hamiltonian system of infinitely many coupled rota- 
tors. We assume that this infinite system is driven out of ther- 
mal equilibrium either because energy is injected by an external 
force (Case I) , or because heat flows between two thermostats 
at different temperatures (Case II). We discuss several possible 
definitions of the entropy production associated with a finite or 
infinite region, or with a partition of the system into a finite 
number of pieces. We show that these definitions satisfy the ex- 
pected bounds in terms of thermostat temperatures and energy 
flow. 
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Introduction. 



In the present paper, we study certain classical Hamiltonian systems consisting of an 
infinite number of coupled degrees of freedom (rotators or "little wheels"). For a system in 
the class considered, the time evolution (/*) is well defined, and given by the limit (in some 
sense) of the Hamiltonian time evolution for finite subsystems. [Note that other infinite 
systems, like gases of interacting particles, would be much more difficult to control]. A 
probability measure on the phase space of the infinite system is called a state, and it has 
a well-defined time evolution. We introduce a family of initial states called F-states (they 
are Gibbs states of some sort). Some of these F-states describe a situation where parts 
of our infinite system (thermostats) are at given temperatures. For a F-state i, the time- 
evolved state gives a finite Gibbs entropy S^{X) to each finite subsystem X of the 
infinite system L. If X is infinite (but has finite interaction with the rest of the system) 
the difference AS\X) = limy^oo('5*(X nY) - S°{X (lY)) still makes sense. 

The bulk of the paper is dedicated to a discussion of the (nontrivial) dynamics of 
our infinite system of rotators. Understanding the dynamics of the system is a necessary 
prereqinsitc to analyzing its nonequilibrium statistical mechanics. We shall in fact ex- 
amine a specific nonequilibrium problem: is it possible to define a local rate of entropy 
production (associated with a finite region X) in a nontrivial manner? This possibility has 
been suggested by Denis Evans and coworkers [16]. We examine their proposal and some 
alternatives, but obtain only partial results. Because of the obvious physical interest of 
the problem, we now give some details. 

By time-averaging f^i or dAS^{X)/dt (over a suitable sequence of intervals [0, T] — > 
oo) we may define a nonequilibrium steady state 



(we do not know that cr{X) is uniquely determined by p and X). 

We ask if an entropy production rate e{X) can be meaningfully associated with a 
finite set X C L. For definiteness we shall think of two physical situations. In Case I 
there is a finite set Xq such that an external force acts on Xq, and the initial state £ 
restricted to L\Xo corresponds to thermal equilibrium at temperature (3~^. In Case II 
we have L = Xq U Li U L2 where Xq is finite, Li and L2 are infinite and £ restricted to 
Li corresponds to thermal equilibrium at temperature /3~^ (with < ^2^)- There is a 
thermodynamic formula for the global rate of entropy production: 




and an average rate of entropy growth 



a(X)= lim -AS^ 



(^) 



cq — (3 X energy flux to thermostat 



(Case I) 



ee = {Pi — P2) x energy flux to thermostat 1 



(Case II) 
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[Note that Case I resembles Case II, where thermostat 2 is replaced by the external force, 
and ascribed an infinite temperature {(32 = 0)]. The question is how to define a local rate 
of entropy production e{X) > such that supj,^^ g^ite ^{^) — ^e- 

The original proposal by Evans and coworkers* is to take, for X finite, 

e(X) = -a{X) 

This is shown to be the average rate of volume contraction in the phase space [X] of the 
subsystem X due to the fiuctuating forces to which it is subjected by the complementary 
subsystem L\X. 

Another idea is to replace the entropy S{X) by the conditional entropy given formally 
by S{X) — S{L) — S{L\X). The corresponding rate of entropy production is 

e(X) = a{L\X) 

We shall make the important physical assumption that the expectation value of the 
energy for each finite system X has a bound independent of time**. It follows that e{X) 
is finite, and one has 

< e{X) < e{X) 

Instead of using a finite set X one may base a definition of entropy production rate 
on a finite partition A = {Xo,Xi, . . . ,Xn) of L, with finite boundary (this will be made 
precise later). We define 

n 

eiA) = J2^iXj) , eiA)= ^ 

j=0 j'-^j infinite 

In particular, in Case II, for X finite D Xq, we have 

e{X) = e{{X, L\X)) < e{{X, Li\X, L2\X)) 

and the right-hand side e{{X, Li\X, L2\X)) seems a rather natural definition of entropy 
production rate. 

We shall later study further properties of the entropy production rates defined above, 
but we note here that they are all bounded by the thermodynamic expression cq. The 

* Actually, the ideas presented in [16] are formulated for Case I, and for a system 
thermostatted at the boundary rather than an actually infinite system. While the two 
idealizations are technically quite different, they are expected to give the same results in 
cases of physical interest. 

** Note that in Case I, if the system has dimension < 2, the external force may cause an 
infinite accumulation of energy in a finite region. Our assumption that the nonequilibrium 
steady state p gives a finite expectation to the energy of finite subsystems is thus invalid, 
and so is our analysis. 
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problem is to prove that they depend effectively on X or A, and are not identically equal 
to or cq. 

We now recall some earlier work to put the problem of defining a local entropy pro- 
duction rate in perspective. 

In earlier studies of quantum spin systems [15], [11], the global entropy production 
(for Case II) was defined by the thermodynamic relation 

ee = {Pi — P2) X energy flux to thermostat 1 

but the quantities e{X), e{X) were not introduced because they would automatically 
vanish. This is because, for quantum spin systems we have < S^{X) (see [3] 

Proposition 6.2.28(b)); for classical rotators by contrast, the entropy is not bounded below. 

The statistical mechanics of classical systems outside of equilibrium can be studied 
in models with nongradient forces and a "deterministic thermostat" [7], [10]. Such a non- 
hamiltonian system corresponds in effect to a rather general time evolution (/*) defined by 
a vector field A" on a finite dimensional manifold M. In general, no absolutely continuous 
invariant measure {i.e., "phase space volume" m) on M is preserved by the time evolution, 
but one may assume that there is a natural (singular) measure p describing a nonequi- 
librium steady state. One can argue that the average phase space volume contraction 
J p{dx){—divmX){x) is the rate of entropy production by the system. This identification 
(for which see Andrei [1]) has been used in particular by Evans, Cohen, and Morriss [6], 
and by Gallavotti and Cohen [9] in the study of fluctuations of the entropy production. 
See also the work of Posch and Hoover [13], Gallavotti [8]. 

Note now that if we introduce a nongradient force ^{q) in the Hamiltonian equations 
of motion, the volume dp dq is preserved, but energy conservation is lost and this is why a 
thermostat is needed. In the case of a deterministic thermostat, the phase space contraction 
is caused by the thermostat (as one can check in the example of the isokinetic thermostat 
corresponding to an added "force" —a{p,q)p, where a{p,q) = p • C{q)/p ' p)- the lab 
however the thermostat is of a different nature: it is typically a large system (reservoir) 
with which the small system of interest can exchange heat, and it is not clear at first how 
to define entropy production. In particular, a nonequilibrium steady state for the infinite 
system L may well have absolutely continuous projection on the phase space of the small 
system X [4], [5], [2], which contradicts e{X) > but may allow e{X) > 0. 

Finally, to indicate the difficulty of the problems considered here, and in particular 
of proving e{X) > 0, consider Case II in dimension < 2. There (as indicated by the 
macroscopic continuous limit), f^i presumably tends to an equilibrium state p and the 
entropy production e{X) vanishes for all X. 

Acknowledgments. 
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correspondence with D. Evans, and many discussions with J.-P. Eckmann, G. Gallavotti, 
J.L. Lebowitz, H. Moriya, H. Posch, and L.-S. Young. 
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1 Description of the model. 

Our system will be an infinite collection of rotators labelled by a; e L, each with 
Hamiltonian H^cipxi Qx) = P^/^ + ^x(fe)) where Px Qx ^ T. [This is for simplicity; it 
would probably be easy to replace the rotators by more complicated systems]. Wc let F 
be a set of unordered pairs {x, y} of points in L, i.e., F is a graph with vertex set L, and 
we define a formal Hamiltonian for the infinite system of little wheels: 

^HxiPx,qx)+ 5Z ^{a;,2/}fe,9y) 

The functions 14, W^{x,y} are assumed to be smooth. 

For X G L, let Fx = {{x, y} eT : x,y & X} and, when X is finite, write 

Hx{px, qx) = HxiPx, Qx) + ^ W[x,y}iqx, Qy) 

where px — {Px)xex G R''^, qx — {qx)xex G T"^. We shall also make use of a constant 
external force* F e R^" acting on a finite set Xq. 

For finite X, a time evolution {fx) on R"^ x T''^ is defined by 

±(px\ fFx-d,^Hx{px,qx)\ 
dt\qx)\ Px J 

where the term Fx is the component of F in R"^, and is present only in case I. We have 
thus 

fx{Px{0),qxm = {px{t),qx{t)) 
We shall suppose that F is connected and, for x,y e L, define 

d{x, y) = min{/c : ^xq, . . . ,Xk & L with xq = x,Xk = y and {xj-i,Xj} e F for j = 1, . . . /c} 

We write then = {y : d{x, y) < k}. 

1.1 Assumption (finite dimensionality). 

There is a polynomial P{k) such that for all x e L and k >0 

\Bl\<P{k) 

[We may take P{k) = 1 + ak^ for some a, 6 > 0; F is thus assumed to have order < a, and 
"dimension" < h]. 

* F is taken constant for simplicity. More generally one could consider the case of a 
smooth function F{qxo, 4>*ct) of qxo and (t)^a with values in R^", where is a smooth 
dynamical system on a compact manifold A, and a is distributed according to some pre- 
scribed (0*)-ergodic measure on A. 
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Note that by compactness of T and the assumed smoothness of T4, W^^-^yy, every 
"force" term {i.e., each component F^; of F for x G Xq, each dq^V^, each dq^W^x,y}j 
dqyW^j:^yj) is bounded, and has bounded derivatives with respect to its arguments q^. 

1.2 Assumption (uniform boundedness). 

The force terms dq^V^, dq.jW{x,y}i ^qv^{x,y} their q^-derivatives (up to any Gnite 
order) are bounded uniformly in x,y G L. 

1.3 Lemma (uniform boundedness of forces). 

The forces F^ — dq^Hx{px-iQx) or —dq^Hxipx-iQx) and their q^-derivatives have 
modulus hounded respectively by constants K, K' . [We shall also denote by K a constant 
> 2K,P{1)K',1]. 

This follows from Assumptions 1.1 and 1.2 [only the bounded order of F is used from 
Assumption 1.1]. [] 

2 Time evolution of infinite systems. 

For X C L, we shall from now on write [X] — (R x T)^. We note the following facts 
which follow from Lemma 1.3. 

(i) For X finite and ^ G [X], if /j^^ = {px{t), Qx{t))xex we have the estimate 

\Px{t) - Px{0)\ < K\t\ when xeX 

(ii) For X finite and f G [X] let also /tf = {Pa:{t), qx{t))xex- Then, ifp^(O) ^Px{0), 
Qx{0) — Qx{0) for some x & X n X we have 

\Px{t) - px{t)\ <2K\t\ 

\qx{t) - qx{t)\ <K\t\'' 
and since \qx{t) — qx{t)\ ^ ^ ^ K, we also have 

\qx{t)-Ut)\<[K\t\^-K]^ <K\t\ 

so that 

max(|p^(t) -px{t)\, \qx{t) - qx{t)\) < K\t\ 
(in) Let /c > and X D B^, X D B^. Then, with the notation of (n), if 

(VyG5^) pyiO)=pyiO) and ^^.(O) = ^^,(0) 

we have 

maxi\p^{t)-p^{t)\,\qS)-Qxm < 
indeed, by the equation of motion and induction on k we have 

\j^Uxi-fxO\<K^^ 
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and the desired result follows by integration. 
2.1 Lemma (a priori estimates). 

For finite X,X C L, let ^ E [X], f G [X], and fj,^ = {Ut))xex = {Px{t),qx{t)), 

fxi={ix{t)),^x = iPx{t),qAt))- With this notation, 

(a) \p,it) - p^iO)\ < K\t\, hence \p,{t)\ < |p,(0)| + K|t| 

(b) ifk > 0, and B^-^ gXUX, and ^y{0) = 4(0) for all y e B^'^ , then 

{K\t\f 



\^x{t) - Ut)\ = maDc[\p^(t)-p^{t)\, \q^{t) - q^{t)\] < 



k\ 



This follows from (i) and (iii) above [(b) is a rather rough estimate, but sufficient for 
our purposes]. [] 

2.2 Proposition (time evolution). 

Let {Px{^),qx{^))xeL e [L] he given. For finite X C L, write {p^ {t), q^ {t))xex = 
fxiPxi^)i 9^(0))x€X- Then for each x E L the limits 

lim p^{t)=px{t) , lim q^{t) = qx{t) 

exist, and (pxit), qx{t))xeL is the unique solution of the infinite system evolution equation 
with initial condition {px{0), qx{Q))xeL- We write {px{t), qx{t))xeL = f{Px{Q), qx{Q))xeL- 

The existence of the limit follows from Lemma 2.1. Writing the infinite system evolu- 
tion equation is left to the reader, as well as checking that ipx{t), qx{t))x€L is the unique 
solution. [] 

2.3 Remarks. 

The limits in Proposition 2.2 are faster than exp{—k d{x, L\X)) for any A; > 0, inde- 
pendently of ipxiO),qxiQ))xeL, and uniformly for t in any compact interval [— T, T]. 

Existence and uniqueness theorems are known in more difficult situations; see for 
instance [12]. 

The proof of Proposition 2.2 does not use the finite dimensionality of F, only its finite 
order. 

2.4 Notation. 

In principle wc use the notation (j)-^ (t), q^ {t))xex for the finite system time evolu- 
tion (fx), and ipxii), qx{t))xeL for the infinite system evolution (/*), but it will often be 
convenient to drop the superscript X. 

It is useful to compactify the momentum space R to a circle R by addition of a point 
at infinity for each x E L, and write [X] = (R x T)^ for X C L. The phase space of 
our infinite system is then [L] = (R x T)-^ C (R x T)-^ = [L]. We shall use the product 



topologies on [L] = (R x T)-'^ and [L] — (R x T)-^; therefore [L] is compact and [L] has 
the topology it inherits as subset of [L] . li U C L we denote by ttu the projection 

7TU : [L] = [if] X [L\U] ^ [if] 

2.5 Proposition (continuity of /*). 

The map t) i— > is continuous [L] x R — > [L] and, for each t, /* : [L] — > [L] is a 
iomeomorpiism. 

To prove the continuity of t) i— > /*^, it suffices to prove the continuity of t) i— > 
{Px{t), Qxit)) for each x G L, and this results from the uniformity of the limits in Proposition 
2.2 (see Remark 2.3). By uniqueness of the map /~* is the inverse of /* and, since 
/~* : [L] [L] is continuous, /* is a homeomorphism. [] 

2.6 Proposition (smoothness of /*). 

Let X CY, X finite and Y finite or = L. For ^ e [X], t] e [Y\X], write (^, t]) = 
{Px{i), (lx{t))xeY ■ Then, for fixed rj and each x & Y, the map {^,t) i— >• {px{t), Qx{t)) is 
smooth [X] X R ^ R X T. 

This results from the bounds on the derivatives (uniform in Y) obtained in Proposition 
2.7 below. 

2.7 Proposition (estimate of derivatives). 
Let Y <Z L,Y finite or = L, and 

fY{Vx{0))xeY = {Vx{t))xeY = {Px{t),qx{t)) 



Write 



where xi, . . . ,Xj need not be all distinct; then 

\ri''^Ht;x^,...,xj)\<Viji{\pym+K\t\)) 

where Vij{{py)) is a polynomial of degree < i -\- 1 (the degree is 1 if (i,j) = (0,0), < i 
otherwise) in the Py with y G B^. . 

Let a = (t(x, xi, . . . , Xj) denote the smallest number of edges of a connected subgraph 
of r having x,xi, . . . ,Xj among its vertices. Then the coefRcients of Vij are positive and 



< MijC 



jK\t\ (-^jl^l 



with suitable Lj, Mij > 0, for all r such that < r < [u — where we have written 
[a — i]+ = max(0, a — i). 

The proof is given in Appendix A.l. 
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2.8 Proposition (estimate of differences). 

We use the notation of Proposition 2.7. Let {r]x{0)), {fjx{0)) G [Y] and define Tx ^ as 
rx'''^with T] replaced by fj. For finite X <Z L we assume r]y{0) = fjy{0) when y ^ X , and 
write 



,Xj;X) = r., 




, . . . ,Xj', X) fx''^^ {t', 



Xi, . . . , Xj'j X^ 



If d{x, X) > i, we have 



Ar^J'^\t;x„... 



,xf,X)\<Q{{\py{0)\ + Km 



where Q{ipy)) is a polynomial of degree < i + 1 (the degree is 1 if — (0,0), < i 
otherwise) in the Py with y E B^. 

Let a — a{x ; X) denote the smallest number of edges of a subgraph of F 

(not necessarily connected) connecting each point x,xi, . . . ,Xj to some point of X. Then 
the coefRcients of Qij are positive and 



for all T such that < r < a — i 

The proof is given in Appendix A. 2. 
2.9 Remarks 

Proposition 2.7, 2.8 will be used in the proof of Theorem 4.5 below. In view of these 
applications the following facts should be noted. 

(a) The condition d{x, X) > i in Proposition 2.8 is not a serious limitation because, for 
the finitely many values of x such that d{x, X) < i, one can estimate Ar*^*'-'^ by Proposition 
2.7 applied to r^^'-?) and f(*'^). 

(b) Write a = a{y,yi, . . . ,yj) and let y be fixed, then 



so that l/{a — i)\ decreases faster than exponentially with respect to ri = \yi — y\, . . . , rj = 
\yj — y\, while \By^ | • • • \By^ \ is polynomially bounded. 



jK\t\ iLj\t\) 
r! 




Indeed, we have \yk — y\ < cr, hence 
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3 Time evolution for probability measures. 

Consider any probability measure I on [L] = (R x T)^ carried by [L] = (R x T)-'^ 
(i.e., £ gives zero measure to the points at infinity). We can find constants > such 
that, if we write 

Bn = {{Px, Qx)xeL : |P:r| < i^nx for all X e L} 

we have (-{Bn) > 1 — 1/n. We may thus write lim^i^oo 11^ ~ = where the measure 
has support in the compact set B^, C [L], and (t, ^ is continuous on R x B^. We 

define then 

fi = hm f£ri (norm limit) 

n — >oo 

Notice also that f^in has support in the compact set B'^ defined like B^ with k^x replaced 
by k'^^ = Knx + K\t\ (see Lemma 2.1(a)). Therefore is again carried by [L]. 

3.1 Proposition (continuity of time evolution). 

If the probability measure £ on [L] is carried by [L] , then the probabihty measure f^£ 
is well dehned, carried by [L], and t i— > f^£ is continuous R — > measures on [L] with the 
vague topology. 

For any continuous function A : [L] — R we have {p£n){A) — £n{A o /*), where 
A o restricted to B^ depends continuously on t with respect to the uniform norm on 
C{Bn — > R). Therefore {f^£n){A) is a continuous function of t, and so is its uniform limit 
t ^ f*£{A). This shows that t i— f^£ is continuous with respect to the w* (=vague) 
topology of measures on [L] , concluding the proof. [] 

Let £x be a probability measure on [X] for finite X G L. We write X oo when, for 
every finite U G L, eventually X ^ U. Suppose that for every finite U and A e C{[U] R) 
the limit 

lim £xiAo7rux) 

exists, where iiux is the projection [X] = [X\U] x \U] \U]. This limit is then of the 
form £{A o tt^) where £ is a uniquely defined probability measure on [L] which we call the 
thermodynamic limit of the £x'- 

£ = e\iuvx^oo^x 

This means that 

'Ku£ = w* limx^ooTTC/X^X 

or (modulo the identifications A Ao tvux, A ^ A o nu) £ is the limit of the £x on 

which is dense in C{[L] — R). In particular, if £ is any probability measure on [L], we have 

£ = 6'limx^oo7rx^ 

We shall later also consider thermodynamic limits associated with a sequence X„ — > oo, 
writing £ = ^limn^oo -^x^ if t^u^ = w* limn^oo 7r[/x„^x„ for all finite U C L. 
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3.2 Proposition (time evolution of thermodynamic limits). 
Suppose that 

e\imx-.oo£x = i 

where £x,^ sue probabihty measures carried by [X], [L] respectively. Then 

dlimx^oofyx = fi 

uniformly for t € [—T, T] . 

We have to prove that, for every finite U C L, and ^ e — > R), 

hm {fxix){A O TTUX) = {fe){A O 7Tu) 

A— »00 

We may (and shall) assume that \A\ < 1. Given e > 0, we know that 

||^o7rj7o/*-^o7r[7xo/^o7rx|| < | 

for sufficiently large X, say X D F for suitable V D U, for all t e [— T, T] . Under these 
conditions we have thus 

\\Ao TTu O — Ao TTUV ° fv ° TTvll < e/2 

and 

1 1^ O T^UX ° fx OTTx - Ao TTjjv O /y O TTyj | < € 

which we shall use below in the form 

\\A o TTux o fj^ -Ao TTuv ° fv° TTvxW < e 
Take now a function $ e C([F] — > R) with compact support and |$| < 1, such that 

\\£ - o 7Tv)i\\ < e 
Using the notation a ~ 6 to mean |a — 6| < e, we have 

{f£){A o TTu) - (/*(($ o 7Tv)mA o nu) = ^(($ o nv){A o nu o f)) 

~ £(($ o nv)iA o TTUV °fv° Try)) = i{{^{A o ttuv o fy)) o yry) = £(*t o ttv) 

where the function ^'^ — ^{Ao ttuv o /y) ^ [^] — * R is continuous with compact support, 
hence extends to a continuous function on [V] . By assumption we have 

|£(*t o Try) - £x(*t o 7rvx)\ < e 
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for sufficiently large X, uniformly with respect to t G [— T, T] (this is because t — > is 
continuous with respect to the uniform norm of C([V^] R-))- We may thus take W D V 
such that, iiX dW aiidte [-T, T], 

e(^t O TTv) ~ ixi'^t O nvx) = ix{{^{A O TTUV o fv)) o T^vx) 

ixii^o7rvx)iAo7rux o fx)) = (fxii^ °T^vx)ix))iAo7rux) 

We have thus 

\{fe){A o TTu) - ifxii^ o 7rvx)ix)){A o ttux)] < 4e 
when X dW, t e [-T, T]. We may now let $ ^ 1, obtaining 

\{f£){Ao7ru)-{fyx){Ao7rux)\<^e 

as announced. [] 

Proposition 3.2 also holds for the thermodynamic limit associated with a sequence 

Xn ^ oo 

4 F-states and their time evolution. 

We introduce now a special set of probability measures. 
4.1 Definition (F-states). 

We say that the probabihty measure i carried by [L] is a T-state if there exist constants 
(3x > (for X e L), smooth functions Vr^ : T ^ H (for x e L) and W[j:^yj : T x T — > R 

(for {x,y} e F; such that the /3^, /S'^ , V^, W^x,y}' ^q.^{x,y}^ 9qy^{x,y} are hounded 
uniformly in x,y & L, and the following holds: 

For every hnite X C L, the conditional measure ixid^lv) of i on [X] given ij G 
is of the form 

£x{d^\v) = const. exp[- ^ {^PxpI + Vx{qx)) - ^{x,y}{lx, Qy)] d^ 

x€X {x,y} 

where J2* extends over those {x, y} & T such that x & X, and we have written ^ = 

{Px,qx)x€X, V = {Px,qx)x€L\X- 

[The F-states are Gibbs states* for a certain interaction given by the Px, Vx, W^x,y}]- 

If £ is a F-state we may, for finite U G L, write {TTu£){d$,) — iuiOdC where £u is 
smooth on [U]. Note that iuiO has a (p, g)-factorization: it is the product of a smooth 

function of the for x G C/, and of a Gaussian sJ^^J^j: exj^(—j3xPl:/'2.) for each x ^ U. 

* See [14] for a discussion of Gibbs states in the simpler case of spin systems. We shall 
not make use of the theory of Gibbs states in the present paper. 
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We shall now take X finite and y = X U Xi U . . . U X^, where 

Xk^{yeL\X:d{y,X) = k} 

is thus the "size" oiY). If ^ e [X], and rjk G [X^] for A; = 1, . . . , fc, the F-state property 
of (. then gives 

■■■,Vk) = CkM^\Vi)-h{Vi\V2) ■ ■■h-iiVk-i^Vk)-hiVk) 
where is a normalization constant and (putting W^^- yj = if {x, y} ^T): 

sex x,y£X xEX y^Xi 

£kiVk\Vk+i) = exp[- X {hvl/'^ + - X] 9y) 

xeXfe x,y€Xk 

X W^{a;,j/}(fe,9j/)] when A;>0, and 

him) = J ^{dVk+i) hiVklVk+i) for some probability measure u on [X^+i]. 
Using the fact that the Jacobian of is 1 (/^ preserves d^) we have 

Thus, by Proposition 3.2, if X C y as above, 

{TTxfim) = w*\imY^oo7rxiiYifY\C,VYmdriY) 

We may write 

fri^, Vy) = ifroi^, Vy), fylii, Vy), • • • , fyH^, Vy)) 

with 

fYoi^,VY)e[X], fyl{^,VY)e[Xk] for A: = l,...,fc 
If I e [X], r]Y e [^\^], the quotient 

£Y{fY'iC,VY)) 

^yUyMvy)) 

is thus a product of quotients 

4(/.-Ke,..)l/.-(W)(e,w)) ^^^^^^ - ^ 
4(/^^(l,w)l/y(Wi)(e~ w)) 
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and 

where the arguments fyli^^VY), /Y(fe+i)(C) their derivatives have dependence on 

^ that decreases faster than exponentially with respect to k (Propositions 2.7 and 2.8). 

Let us define '^ykiC-iVY) by 

4(/y^(e,w)i/y(Wi)(e,w)) = 4/c(e,w)-exp y1 {-kp.i^fm 

for k — 0, . . . ,k — 1, where X}^ is replaced by X for A; = 0, and 
We shall also use rj) defined by 

xeXk 

where 

with fo\^,v) e [X], and fj;\^,v) e [Xk] for /c > 1. 
From our definitions it follows that 

4.2 Lemma (basic uniform estimates). 

In the above formula, we have, uniformly in t e [— T, T] and the size k of Y , the 
estimates 

(a) I log4o(C: Vy)\ < const. (1 + sup \px{0)\) 



W |logf«|<^^(l+ sup \M0)\) if^>l 

These estimates remain true when ly^ is replaced by l\. 
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We note that, by Lemma 2.1(a), 

\Px{tf-Px{Of\ < K\t\{2\p^{0)\ + K\t\) 

From this, and the definitions, the first inequahty of the lemma follows. The second 
inequality is obtained by using also the finite dimensionality Assumption 1.1 and Lemma 
2.1(b). D 

Define now the regions Ru, Ry C [L] = [X] x such that 

Ru = {it T]) : \p^\ < u if X e X} , R^ = {(^, ry) : \p^\ < kv ii x e X^ for k > 1} 

4.3 Lemma (existence of limit in R^ n Ry). 
In Ru n Ry , the expression 

r ^^Ji\f-'f ^^ " [4o(e, w) exp ^ (-^(0)72)]"^ 



/ "m 

J\X] 



I "-11^* (, ..4o(l,w)exp^(-^,p,(0)V2) 

'[X] i^^.^YkitVY) ^x 

has upper and lower bounds cxp(±const.(l + v)) uniformly in u, t & [—T,T], and k, and 
tends when k —>■ oo, uniformly for (^, ry) € Ru H R^ , to 



?o(e,^)exp5^(-4px(0)V2) 



'[X] t^.m^v) 

The limit is continuous. 

Let (^, 77) e -Ru n Ry , and assume k to be large. The quotients 

^YkilvY) 



'Yk 



itVY) 



ik>l) 



arc nearly independent of Y (i.e., of k) for small k, and (using Lemma 4.2) very close to 
1 for large k, so that 



-Yk 



(I, Vy) 



n 



uniformly, and we have bounds exp(±const.(l + v)) by Lemma 4.2(b). Note now that 
■^yo(C) Vy) tends to ^0(^5 v) uniformly for (^, 77) e i2„ fl Ry , and we can extend the integral 
over ^ from \px\ < u to [X] because the Gaussian 

exp J](-^,p,(0)V2) 
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beats the exponential growth of ^yq{^,i]y) given by Lemma 4.2(a). Bounds of the form 
exp(±const.(l + v)) hold again after integration. [] 

4.4 Lemma (large v Gaussian estimate). 

For large v, (7ry\x/r7ry£)((i?7y) has mass < exp{— const. v'^) outside of tiy\xRv ' 
formly in the size kofY. 

The (p, g)-factorization of £y{^-,Vy) shows that the mass outside of is bounded, 
uniformly in k, by a Gaussian < exp(— const. v^) for large v. But the time evolution fy 
changes \px\ (additively) by at most K\t\, so that the Gaussian estimate remains valid. [] 

4.5 Theorem (Smooth density of evolved states). 

Let £ be a T-state. For finite X, and Y of size k as above, we write 

^lx(Oexp5^(-te(0)V2)= / driYiY{fYi^,VY)) 

There is a smooth function ixiO of C and t such that 

(nxfim) = 4(6 exp J2 (-te(0)V2) 

x€X 

and we have, uniformly for \px\ < u {x E X) and \t\ <T, 

4(0 = _lim £1^(0 

fc— »oo 

The limit also holds for the derivatives with respect to ^, t. The /^(C)? ^yx(0 have upper 
and lower bounds exp(±const.(l + u)), and the absolute values of their derivatives have 
bounds polyn.(tt). exp(const.(l + u)) uniformly in t e [— T, T] and k. 

We start with the remark that 

^Y{fY\C,VY))d^ 



Ilx]d^MfY\i^VY)) 



is the conditional measure of /y'^y'^ on [X] given r]Y G . Integrating this conditional 

measure with respect to {7rY\xfY'^Y'^){dr]Y) yields ttx/y'^y'^- Thus 

^lx(e) = / {TTYXxfl^TTYeXdVY) , ^lli^'^^t^P " exp J] (te(0)V2) 

J I[X]d^£Y{fY\^,VY)) ^x 

The integrand in the right-hand side is the product of a factor controlled by Lemma 4.3, 
and a factor ^yolC^ 'Ty) which has upper and lower bounds exp(±const.(l + u)) uniformly 
in k (by Lemma 4.2(a)) and tends to £o{^, rj) when ^ — > oo, uniformly for (^, r/) e i?^ fl . 
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Using also Lemma 4.4 and the fact that TTyxxfy'^Y'^ has the t(;*limit 'KL\xf^^ when A; ^ oo 
(Proposition 3.2) we find that 

Jim £1^(0 



= {^L\xtmv)m.v) 



°° ot 



uniformly when |px(0)| < w for x G X, with uniform upper and lower bounds exp(±const. 
(1 + u)). We call the limit lx(0- Since 

^Yxii) exp 5] (-to(0)V2) = (TTx/f 7rr^)(rfe) 

has the w* limit {7rxf*i){d^), it follows that this limit has a density 

4(Oexp5^(-te(0)V2) 

as asserted. 

Using the notation 



we may write 



^^~[n li4]^o(l,^)exp 5^(-4p.(0)V2)l"'^^(e,^) 



[X] ' ,iVfc^^'^) 



xex 



= / ^(c^^cZry) 



[H iii^j\irj)/ '^^^ ^K^^^ ^ (-4p.(0)V2) 



xex 



and remember that this is the limit of a similar expression for ^xy(0- want to 

show that derivatives (of all orders) with respect to $,,t by showing that the 

derivatives of i^xviOi in a compact set, are bounded with respect to k. Note that 

in estimating the integrals of polynomials in p, the Px{0) integral always has a Gaussian 
factor exp{— PxPj:{0)'^ / 2) (remember the (p, q) factorization of £). Therefore we only have 
to worry about bounding the coefficients of the polynomials. 

Inspection of the above expression shows that computing a first order derivative es- 
sentially involves multiplying the integrand by the logarithmic derivatives of /o(^, f^i'n))^ 

or ^fc(^, /|(?7))/^fc(^, /|(?7)) and summing over k. In view of the very explicit form of the 
logarithm of the £k, we just have to estimate the derivatives of /~*(^, /|(^)) with respect 
to ^ and t. 

As far as ^o(^, f^{f])) is concerned, we have to consider the derivatives of /o~*(^, /|(^)) 
(or /i~*(^, /|(^7)), which is similar). The ^-derivative is of the form r^'^\—t;xi) with 
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x,xi e X, giving a bounded contribution by Proposition 2.7. The t-derivative is of the 
form 

with X E X, and we may take Xi G Xk-^- By Proposition 2.7, |r*-°'^-'(±t)| is bounded by a 
polynomial in p with bounded coefficients, and r^'^\—t;xi) < const. //ci!, again giving a 
bounded contribution because J2ki l^fej/^i' is bounded. 

We turn now to -^felC) /|(^))' '^^ have to consider the derivatives of 

/fc"*(C, /|(^)) (or The ^-derivative is of the form r^x'^\-t; xi) with x e Xk, xi e X, 

so that \r^'^\—t;xi)\ < co-ast./d{x,xi)\ = const. /A;!, which makes a bounded contribution 
because \^k\/k\ is bounded. The t-derivative is of the form 

Ar(i'°) {-t; X) + J2 Ar(«'^) {-t; x,; X)ri\''^ (t) 

Xl 

where x e X^ and we may take xi e X/-^. By Proposition 2.8 we have 

\Ar^J'^\-t;X)\ <polyn.(p)/A;! 

\Ar^^'^\-t;xi;X)\ < const. /o"! < const./ max(A;, fci)! 

and since \Xk\/k\, l-^feM-^fei 1/ max(/c,A;i)! are bounded, we also have bounded 

contributions for the t-derivative. 

We consider now higher order derivatives with respect to ^,t. The computation of 
such a derivative gives terms where the integrand is multiplied by a product of logarith- 
mic derivatives of the type discussed above; the contribution is again seen to be bounded. 
There are also terms containing derivatives of the logarithmic derivatives, and these are 
expressed in terms of higher order derivatives of with respect to ^,t. The 

derivative /d^xi • ■ • d^xj is estimated by \r^'^\—t; xi, . . . , Xj)\ < const. /d{x, X)\ giving a 

bounded contribution. The derivative /df is a sum of terms Arx^'^^\—t; yi, . . . , yi^;X) 
(multiplied by derivatives fydt^ = r*^^''^^), where i\ -\- 12 = i', these terms can be esti- 
mated by Proposition 2.8, and give a bounded contribution. The general mixed derivative 
d^+^ /dVdix^ - ■ -di^., with j > 1, is a sum of terms ri*''^"^*^''(-t; xi, . . . , x^-, yi, . . . , y^J 
with xi, . . . ,Xj E X (multiplied by derivatives of the form r'^'^'^^) which can be estimated 
by Proposition 2.7, and give a bounded contribution. [] 

4.6 Remark (uniform bounds). 

The proof of Theorem 4.5 gives estimates of ^yxi'O ^^id its derivatives with respect 
to t and ^, which are uniform with respect to the size A; of y. They are also uniform with 
respect to the P-states £ with conditional measures corresponding to a fixed choice of 
Vx, W^x,y}j cind remain uniform if some of the W^x,y} ^ire replaced by 0. 
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For Y finite and rj = (px, qx)xeY define 



■2' 

xeY x,yeY 



where Z^^ is a normafization factor, and write 



xex 

Then the above remarks show that the uniform estimates on its derivatives 

given by Theorem 4.5 can be taken to hold also for lyx- 

5 Entropy. 

Given a F-state £, and X cY finite, we write 

where {t, ^) i-^ ^YxiO^^^xiO smooth on R x [X]. 
In Theorem 4.5, we used the notation 



^YX 



(0 = ^YxiO exp J2 i-M^fm , 4(0 = 4(0 exp (-te(0)V2) 

xex xex 



and saw that iyxiO tends to •^x(0 together with its derivatives, uniformly on compacts 
of R X [X], when Y ^ oo. 

We can now define a (Gibbs) entropy Sy{X) or S*{X) by 

s'y{x) = - / 4^(0 iog4x(0 , s\x) = - f 4(0 iog4(0 d^ 

J[X] J[X] 

These are convergent integrals in view of the uniform bounds given in Theorem 4.5. Fur- 
thermore Syi^X) S*{X)j uniformly for |t| < T, when Y oo. 

We may assume that Y D Xi ^ {y E L : dist(|/,X) = 1}. If ^ G [X],ri G [Y\X] 
or let rji e Xi be obtained from rj by restricting the index set to Xi. Then the 

equations of motion for ^, 77 show that we may write 

where X does not depend on Y. Writing £^ = iy o we have 

4x(0= / dfjiU^r^) 

J[Y\X] 
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and the "continuity equation" 



It 



so that 



log 4^(0 = -TT^V^ ■ / dr^eUC,rJ)X{^,rJ,) 

:^[4x(e)iog4x(e)] = -(iog4x(e) + i)v^ • / dveuc,v)xic,vi) 

J[Y\X] 

Using the estimates of Theorem 4.5 we find that t SyiX) is a smooth function of 
t, with 



d 
dt 



SUX)= deiog4x(0Ve- / dr7£l(e,r7)A'(e,r7i) 
J[x] J[y\x] 

J[x] %xUj 

Write now X+ — X U Xi = {y e L : d{y, X) < 1}. The probability measure 

^Yx+i^,Vi)dCdrii 
conditioned on ^ e [X] is denoted by iyx+iVilO diji where 

Theorem 4.5 gives uniform estimates for 

YX+ 

so that 

J[x] i[Xi] 

= / c?e4x(0Vr / c^^i4x+(^iiO'^(e,r/i) 

It follows also that, when Y —>■ oo, dSy^X) / dt tends to 

[X] J[X^] 
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uniformly for \t\ < T, and the limit is dS*{X)/dt. 

5.1 Proposition (time derivative of S{X), X finite). 

When Y ^ oo, the derivative dSy^X)/ dt tends, uniformly for \t\ <T, to 



d_ 
di 



J[x] J[x^] 



which is a smooth function of t. 

The proof, as given above, is essentially a corollary of Theorem 4.5. [] 

Suppose now that Xi = {y e L : d{x, y) = 1} is finite, but X is not necessarily finite. 
We still have, for Y finite. 



lot 

dt 



sUxnY)^ / dc x+nY 

J[xnY] J[Xi] 



'[xnY] J[Xi 

and this can be bounded independently of Y. 

5.2 Proposition (time derivative of AS^{X)). 

IfXcL, and X is not necessarily finite, but Xi = {y e L : d{X, y) = 1} is finite, we 
may define 

AS\X) = lim {S\X nY) - S^{XnY)) 

Y— >oo 



and we have 



^AS'y{x)= f [irxfmovr I dm£'x+imm{^,m) 

J\X] J\XA 



'[X] J[x^ 

which is a smooth function of t. 

This follows from the usual estimates. Note that is a derivative with respect to a 
finite number of variables corresponding to nonzero components of X{^,r]i). [] 

We shall now study a conditional, or "external" entropy S defined for X finite by 

S'y{X) ^ S'yiY) - S'y{Y\X) 

Using the notation 

i'Y{^,v) 



4 = iY o fy' , iY,Y\xiv) = / d^^Ut V) , ^Yxi^lv) 

J[X] ^Y,Y\X 

as above, we find 



iv) 



tr. 4 



SUX) = - / d^dvfy{^,v)log 



^Y,Y\X{V) 
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where we have written 



S'y{X\^) = - J d^fyxm^Oge^y^ilv) 



Let also 

J[x] f^^^^kK^^V) 

x£*(e,r7)exp^(-te(0)V2) 
xex 

S\X\rj) = - [ d^i'xm^ogfA^ 



then Lemma 4.2(a) and Lemma 4.3 imply that Sy^Xlr]) S^{X\r]) when Y — > oo, uni- 
formly for 7] G TtL\xRy and t e [— T, T], and with uniform bounds 

\S^{X\r])\ < const.(l + 'i;) 

Therefore (using Lemma 4.4 and Proposition 3.2)we see that when y — > oo we have 

S^y{X) S\X) 

uniformly for t e [— T, T] , where 

S\X) = J{7rL\xfi)idv)S\X\r,) 

Note that <S'*(X) is obtained by taking the mean entropy S''^{X\r]) associated with i^xi^lv) 
and averaging over rj, while S^{X) is the entropy associated with the average •^x(C) of 
•^x(^l^) ^^^^ V- III particular, concavity gives S^{X) < S^{X). 

Since dSY(Y)/dt = 0, we have 



[Y\X] J[X] 



dr] 



^*Y,Y\x{v) T^^vy- [ d^^\.xmyy{^,v) 

~. J\x\ 



^^^^ .[X] 

where X\ = {y : d{y,X) = 1} and yy is the j/-component of 3^. We may now let 
y — > oo, finding: 
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5.3 Proposition (time derivative of the entropy S). 

When Y ^ oo, the derivative dSY{X)/dt tends, uniformly for t e [— T, T], to 



dt 



vex, -t^l 



which is a smooth function of t. 

This is again a corollary of Theorem 4.5. [] 
Note that 

5.4 Assumption (bounded energy). 

For every finite X <Z L the kinetic energy is bounded independently of t: 



c^e4(0Px/2<const.(X) 

[it would be equivalent to assume a bound on the total energy Hx] 
We have the general inequality 



S\X)< j diex{i)p\/2 + \X\\ogV2^ 



[this follows from the "variational principle for the free energy", and can be proved by 
using the concavity of the log: 

j rfC4(0log-^^ < log j d^e-P-/^ = |X|log^ ] 

Therefore the bounded energy assumption gives a bound on the entropy: 

S\X) < const. (X) 

Similarly, we find 



hence 



S\X\r^)< y"de4(eWPx/2 + |X|logx/2^ 
S\X) < j dCixiOPxf^ + \^\ log < const. (X) 



In particular we have S\X) < S\X) < const. (X). 
5.5 Definitions (large volume limit). 
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1 f " 

We may take a sequence (T^) tending to +00 such that — / dt /*£ has a hmit p 

Tn Jo 

in the vague topology of measures on [L] : 

1 Z"^" 

J-n Jo 

We call the probability measure p a nonequilihrium steady state (NESS). In view of As- 
sumption 5.4, p is carried by [L\. Furthermore p is invariant under (/*). 

We can also (by going to a subsequence) assume that 
A5'^"(X) 1 , d 



Tn 



when Xi = {y E L : d{X,y) = 1} is finite (X need not be finite). Note that cr(X) might 
not be determined by p and X. 

For notational simplicity we shall write T — > 00 instead of T„,n —> 00. 

5.6 Interpretation (entropy production). [16] 

As mentioned in the Introduction, Denis Evans and coworkers [16] have proposed to 
identify the mean entropy production rate in a finite region X to 

e(X)=-.(X) = -lin.^"'^)-^°W 

According to Proposition 5.1 this is the mean rate of volume contaction in [X], and e{X) 
corresponds to the accepted definition of entropy production in the presence of a deter- 
ministic thermostat. 

A related choice is 

e(X) = a(L\X) = - ^lim ^ 

This is the mean rate of volume expansion in [i^\X] , and corresponds to the rate of entropy 
growth due to X, as seen by the "external world" L\X. 

Since S^X) < S\X) < const. (X), we have < e{X) < e{X). 

We may also define mean entropy production rates associated with a finite partition 
A = {Xq, Xi, . . . ,Xn) of L provided Xo,Xi,...,X„ have finite "boundaries" {y e L : 
d{Xi,y) = 1}. We write 

n 

e{A)=J2^{X,) , eiA)= ^^^3) 

j=0 j'-^j infinite 
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In particular, in Case II, for X finite D Xq, we have 

e{X) = e{{X, L\X)) < e{{X, Li\X, L2\X)) 

We proceed now with some general inequalities satisfied by cr, e, and e. 
5.7 Basic inequalities. 

We have e(0) = e(0) = by definition, and remember that < e{X) < e{X). The 
strong subadditivity of the entropy implies that, if U, V have finite boundaries, 

a(U UV)- a{U) - a{V) + a{UnV)<0 

(we have used the fact that S^{{U UV)nY) - S^{U nY) - S^{V nV) + S^{U nV f] 
Y) is bounded independently of Y). This implies the strong superadditivity of e, and 
subadditivity of e. In particular 

e{UUV)>e(U) + e(V) if UnV = $ 

and 

e{UUV) < e{U)+e{V) 

If [/ C F we also have* 

eiU) < e(V) , e(U) < e{V) 

i.e., e{X), e{X) are increasing functions of X . 

We can extend the definition of e{X),e{X) to infinite X: 

e{X)= sup e{U) , e{X) = sup e{U) 

finite UcX finite UCX 

In the situations of interest for us e{L) will be finite and we may call this quantity the total 
entropy production rate. Note that the entropy production rate e{X) is not an additive 
function of X, but that its subadditivity amounts to some kind of locality. Note also that 
if e{X) > we must have S^{X) —oo, in particular the ij^ cannot remain bounded 
when T — > oo, contrary to some evidence [4], [5]. But there is no obvious objection to 
having an entropy production rate e{X) > 0. 

6 Thermodynamic bound on entropy production. 

We shall show that in Case I (an external force and a thermostat at temperature 
we have 

e{X) < P X energy fiux to thermostat 

* Note that, for U GV,we have S^{Y\U) < S^{Y\V) + S^{V\U), hence S^{V) < 
S^{U) + S^{V\U), hence e{V) > e{U) + e{V\U). 
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where the right-hand side is the thermodynamic rate of entropy production. A more 
general result is given below (see Proposition 6.3). 

In Case I we have introduced a finite set Xq on which external forces act. As initial 
state £ we shall use the thermodynamic limit of a sequence: 

£ = 9 limy^oo^y iv) dv 

where* 

£y{v) = Zy^ exp[-#Xo - PHy\Xo] 

In this formula, 

xeY\X x,yeY\X 

and Zy^ is a normalization factor. In this section it will be convenient to use Xq instead 
of X in the definition of Y, so that Y = Xq U Xi U . . . U X]^. We take X of the form 
Xq U Xi U . . . U Xk_ (this is no serious restriction) and choose a subsequence A; — > oo 
such that the 7rxY{£Y{'n) drj) converge vaguely (we use here the thermodynamic limit for 
a sequence as explained in Section 3). 

The state £ is a F-state corresponding to the choice Px = = f3Vx for x ^ Xq^ 

and = PW{x,y} for x^y ^ Xq. We define 

4x(e)= / ^Y{fY\i,rjY))drjY 

J[Y\X] 

and 

~s'y{x) = - [ 4^(0 log £1^(0 rfe , Kix) = sUy) - sUy\x) 

J[X] 

Writing 

^YXm=^YY{^,V)/hY\x{v) 

S'y{X\v) = - j dftyxm\0g~eyx{i\r^) 

we find 

Sy{X) = [ £ys^Av)dvS'yiX\v) 



* More generally we could allow a term ^^.^Xq y^Xo ^{a^.J/} interaction between Xq 
and Y\Xq. 

26 



6.1 Lemma (thermodynamic limit for S). 



together with the t-derivatives, uniformly for t e [— T, T], when k —>■ oo. 

We have shown (in the proof of Proposition 5.2) how Syi^X) — > S^{X). We proceed 
in the same way here, using the uniform estimates of Theorem 4.5 which hold again when 
£ is replaced by £, as explained in Remark 4.6. [] 

Since fy is volume preserving in \Y] we have 
therefore 

sl{x) - Sy{x) = sUy\x) - sUy\x) 

We fix now X, with Xq C X cY as indicated above. Note that, by the (p, q) factorization, 
Sy{Y\Xo) = Sy\Xo{Y\Xo) is the sum of a momentum term S^p (integral over p, trivial) 
and a configuration term S^'^ (integral over q). The configuration part yy;5^(9y\x) of 
^Y,Y\x{qY\x) differs from £y,y\Xo{(1y) = £y\Xo(^x\^oi(lY\x) by a factor bounded inde- 
pendently of Y (because there is a finite number of bounded terms and W^x,y} with 
X e X\Xo). Therefore | log£yY\x(^y\x) ~^^S'^y\Xo^^^^\ bounded independently of y, 
hence \S^{Y\Xo) - S^{Y\X)\ < Cq with Co independent of Y. Define now a function £* 
on [Y\Xo] = [X\Xo] X [Y\X] to be the product of 

Z-'eM-PPxxxof^) 

on [X\Xo] (with a normalization factor) and iyy^x [^\"^]- Then there is a 
constant Ci such that 

S* = - [ drj£*{ri)\og£*{ri) = S'y{Y\X) + Ci 

J[Y\Xo] 

and the "variational principle for the free energy" gives 

S* - SUy\Xo) < f drj (r (r/) - £y\xM) PHy\Xoiv) 

SO that 

Si.{Y\X) - S^{Y\X) <Co-Ci + S* - S^{Y\X) 
<Co-Ci+ [ dry (r (ry) - £y\Xo{v)) P Hy \Xoiv) 

J[Y\Xo] 
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There are also constants C2, C3 such that 

j[Y\xo] J[y\x] 

f dri£Y\Xo{v))PHy\Xoiv)> [ dviy^x{v)/3Hy\x{v) + Cs 

J[Y\xo] J[y\x] 

Therefore, with a constant C = Cq — Ci + C2 — C3 independent of Y and t, we have 



sUx) -sUx)<C+ f dii[£'y^y^Av) - ^'Y,Y\xiv)my^x{v) 

J[Y\X] 

= C + p [ drj [4(r/) - iUr})]HY\x{r}) 

J[Y] 

= C + p[ drriUv)[HY\x{fi^v)-HY\x{v)] 

J[Y] 

= C + pJ^ dviUv) jyr-^HY\x{f^v) 



The equations of motion yield 

A 

d^^ 



HY\x{fYV) = ^{7^X+,YfYV) 



where X+ = {x & L : d{x, X) < 1} and , gy) is given by 

d 

xex y^x 



Therefore 



/ dry £1(77) f dT^7Tx+,Yf^v)= f dr [ drjtY{vM7Tx+,YV) 

J[Y] Jo Jo J[Y] 

= fdr f d^iW+{i)m 
Jo J[x-^] 

so that ^ 

Sy{X) -Sy{X)<C + P f dr f df^Yx+mi^) 

Jo J[x+] 

We may now let F — > 00, obtaining 

s\x)-s\x)<c + (3 fdr [ d^e^x+iOm 

Jo Jlx+] 
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6.2 Proposition (bound on entropy production, Case I). 

In case I the mean rate of entropy production of the finite set X is < f3xenergy Rux 
out of Xq: 

< e{X) < p j i7rx,+pm)M0 

where $o is the function $ computed for X — Xq. 

It suffices to consider tfie case of large X, so we assume X D Xq. Taking in the 
previous inequality the large time limit described in 5.4 we obtain 

0<e{X)<P J{nx+pm)Hx) 

where the right-hand side is independent of X, and we may thus take X — Xq. [] 

We now give without proof a general bound on (t{X), which can be obtained using 
the same ideas as for Proposition 6.2. 

6.3 Proposition, (bound on a{X)). 

Let X be infinite, with finite "boundary" Xi = {y & L : d{X,y) = 1}. we let the 
initial state £ be the thermodynamic limit of a sequence: 

£ = ^limy^oo^y (??) drj 

where 

irirj) = Z-^e-^- 

X&Y x,y£Y 

We assume that 

Px = P , V, = PV, , = 
when x,y E X , i.e., X is a thermostat at temperature j3~^ . Then 

(t{X) < P X energy flux to X 

Note that if the energy flows out of X, then a{X) < 0, and in particular a{X) does 
not vanish. Applications of Proposition 6.3, in particular to Case II, are left to the reader. 
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A Appendices. 

A.l (proof of Proposition 2.7). 

By uniformity of the bounds it suffices to consider the situation when Y is finite. The 
case = (0, 0) follows from Lemma 2.1(a), with r = a = 0. 

Differentiating the time evolution equation for ^^^{t) = ipx{t), qx{t)), we find 

where ^xy depends smoothly on the g's, is independent of the p's, and there is a uniform 
bound 

y 

Therefore, if r{t) = sup^, ^.^^y \r^'^\t;xi)\, we have 

||r(t)|<Kr(t) 

with r(0) = 1, so that \r^x'^\t;xi)\ < r{t) < e-^l*l. 
With the notation k — d{x, xi) we claim that 

\ri'''\t;x,)\<{e^\% 

where we have defined {e^}k = — X]n=o ^"/^' ^ inf^<A; u^/£l For k = the claim has 
been proved above. For /c > we have by induction 

\j/^''Ht-x,)\<K{e^W},_, 

and our claim follows by integration, using r^'^\o, xi) = 0. 

Remember that a{x, the smallest length of a subgraph of F connecting 

We claim that for j >1 there are constants Lj > and Mj > such that 

|r^)(t;xi,...,x,)|<M,e^-^l*l inf (1) 

0<T<a T\ 

where cr = (t{x, a;i, . . . , Xj). 

We have already proved (1) for j = 1, with Li — K, Mi = 1. For j > we have 

-^r(°'^)(t;xi,...,x,)= ^ ^xyri''^^\t;xi,...,Xj)+rest (2) 



dt 
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The rest is a sum, over y e and X, of products of factors r^Jl'^"'\t; Xn) where Zn is x 
or y, X = (Xn) is a partition of (xi, . . . , Xj) into \X\ > 1 subsequences of length jn, and 
each product has a coefficient which is a smooth function of qx,Qy Thus 



\rest\ < Cjmaxexp(^j,K|t|)J] ^^^"'y"" < Cjc^^^-^ML 

n n 

where L'^ = max^r ^„ Lj„ and r' must be of the form Tn, i.e., < r' < o-(zn, Xn) 
where each Zn is either x or y, and therefore 

(^(^^n, ^n) + 1 > Cr(a;, Xi, . . . , X^) = CT 

n 

SO that aU values of r' between and [u — 1]+ are allowed. We have thus 



^ 0<T'<[a-l]+ t'\ 

We shall prove (1) by induction on j, assuming now j > 1- First let us write 
sup \r^^'^^ {t; xi,...,Xj)\^r{t:xi,..., xj) = e^^*^s{t) 

X 

and let 

a' = cr{xi, . . . , Xj) = mincr(a;, xi, . . . , Xj) 
Then, for < r' < [cr' - 1]+ and t > 0, 



— I 
dt 



-r{t; xi,...,Xj) < Kr{t; xi,...,Xj) + Cje^^^^— 

T ! 



or 



dt ' ' - ' t'\ 
Thus 

^ < d^e^J-mt^m^^ for 0<r'<a'-l 
dt - dt^ (r' + 1)! ^ ~ ~ 

and also (taking t' — 0) 

ds ^ d Cj ^ ^(j-i)Kt 
dt - dt (j - 1)K 

We shaU take Lj = L'^ + Cj and Mj = 1 + Cj/L'j + Cj[{j - 1)^]"^ In particular, we have 

sit) < M.e^^-^^^*^-4—— for < r' < a' - 1 
\ ^ - J (r' + 1)! ~ ~ 
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s{t) < Mje^^-^^^' 

hence 

< M,e(^-i)^* inf 

o<T<a' r! 

and 

0<T<a' T! 

To prove (1), it suffices to show that if 

a{x, xi, . . . , Xj) > a' + k 

then 

iLj\t\y 



\r^^'^\t; xi, . . . , Xj)\ < M^e^^\^\ inf 



■' ^ o<r<(T'+fe r! 

and wc have just shown this for A; = 0. We proceed now by induction on k for /c > 0. For 
y & we have 

a{y, xi, . . . , Xj) > a{x, xi, . . . , Xj) — l>a' + k — 1 
Therefore, our induction asssumption and estimate of \rest\ give 

dr^^'^){t;xi, . . .,xj)\ < KM,e^^\'\ inf il^ML. + c.e^m i^f (^ML 
Since Lj = + Cj > Cj + and 1 < Mj, we may write for t > 0, 



hence 



, \r^^"-\t;xi, . . .,xj)\ < Mje^^\'\ inf i^,-^ 



|r^•)(^;a;l,...,a;,)|<M,e^•^I^I^M: 

if 1 < T < cr' + /c, but the above inequahty also holds by the induction assumption when 
T = 0, and this completes the proof of (1). 

We discuss now the case i > 0. We have an explicit expression 
given by the evolution equation for {p, q) if j — 0, by (2) if j > 1. We may differentiate 
repeatedly with respect to t, replacing the derivatives in the right-hand side by using either 
the evolution equation for (p, or (2). We express thus r^*'-'^ as a polynomial in the py 

and the ry^'^"* with I < j, with coefficients that are smooth functions of the qy. The indices 

y of the Py,qy, and r^'^^ that occur satisfy d{x,y) < i. Furthermore, in any given term 

of the polynomial, the factors rijl'^^\Xi), . 
forming a partition of {xi, . . . ,Xj)) satisfy 

a{x,yi,...,ym) < i 
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of the polynomial, the factors rijl'^^\Xi), . . . , rylf"^\Xm) that occur (with {Xi, . . . , X^) 



a{x,yi,...,y„i) + a{yi,Xi) + . . . + a{ym,Xm) > (t{x,xi, ...,Xj) = a 

so that 

+ • • • + a{ym,Xm) >a-i 

Therefore, using (1), we see that r^'^'^^ is a polynomial of degree < i in the Py such that 
y E Bl., with coefficients bounded in absolute value by 

const.e^^^l*! inf ^J^ML 

0<T<a-i r! 

if z < cr, otherwise by 

const.e^-^l*! 
concluding the proof of the proposition. [] 
A. 2 (proof of Proposition 2.8). 

Note that the conditions on the coefficients of Q are of the same form as those on the 
coefficients of V in Proposition 2.7, but a has a new definition (and the L^, may have 
to be chosen larger than for Proposition 2.7). 

Note also that Lemma 2.1(b) gives, for d{x, X) > 0, 

|AriO'0)(t;X)| < if l<T<d{x,X) 

Using also Lemma 2.1(a), this proves the Proposition in the case = (0,0) since here 
cr = d{x, X). 

We shall later use the fact that for the qf-component we have actually, for d{x, X) > 0, 

(K\tW 

\Aqa,it;X)\ = \q,it) - ^,(t)| < if < r < dix.X) 

r! 

[this is because qx{t)iQx{t) G T, so that \qx{t) — qx{t)\ < 1]. 

In the study of Ari°'"''* for j > we do not impose the condition d{x^ X) > 0. If j > 0, 
we claim that there are constants Lj > and Mj > such that 

\Ar^^'^\t;x„...,xj;X)\<Mje^^\\M (3) 

0<T<<7 t! 

where cr = cr(a;, Xi, . . . , Xj;X). 

This will be proved by induction on j, using (1) and the equation 

^ Ari°'^)(t; xu...,Xj;X)= ^ $,j,(g„ qy)Ar^'>'^\t; . . . , x,; X) + rest 



dt 
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The rest here is a finite sum of products, each of which has exactly one factor with a 
A in front of it. The factors are: a coefficient depending smoothly on qx,qy, and factors 
r^^'-'"-' where is x or y and the form a partition X = {X^) of {xi, . . . ,Xj) into \X\ 
subsequences of length 

In particular, for j = 1, the rest is A^^y-fy''^^ ■ Using the remark above on [Agfa;!, 
and the earlier bound on jri^'"'^''! one finds 

(K\t\)'' KifiiK\t\y ^|,|(2K|t|)*+^ 
rest] < const. ^4r^-e ' 'Mr^ < const.e^'^l \, ' ' 

where k is allowed values in [0, d{y, X)] or [0, X)] and £ is allowed values in [0, d{y, xi)], 
so that k + £ is allowed all values in [0, [a — !]+]■ 

For general j > 0, using induction on j, and the bounds on jAg^:!, shows that 

the products appearing in the rest have, in absolute value, bounds of the form 

_,(S|)!.(expi:,„A-w).n^ 

n n 
ik + J2£ny. - T'l 

where Lj = K + max^;- Lj^, and we must now discuss the range of r' = k + Y^£n- 
Remember that there is a A in front of one of the factors of the product we are considering. 

If the A is in front of the coefficient depending smoothly on qx,qy, this corresponds 
to k E [0, d{z, X)] with z = x or y, while in G [0, (T(^n, -^n)] with Zn = x or y. Since 
d{x,y) = 1, we have d{z,X) + J^i'^i^e^ ^e) + 1 > a{x,xi, . . . ,Xj;X) = a; therefore 
t' = k + ^ £n is allowed all values such that < r' < [o" — 1]_|_. 

If the A is in front of one of the ri°'"'\ say for n = a, the corresponding is G 
[0,a{za,Xa:,X)] by the induction assumption, the other ri^'"'"'* are e [0,a{zn,Xn)], and 
we have A; = 0. Note that 

a{za, Xa, X) + ^ a{zn, Xn) + 1 > a{x , Xi , . . . , Xj , X) = a 

therefore r' = is again allowed all values such that < r' < [u — 1]+. In conclusion 
we have 

\rest\ < Cje^^\'\ inf ^^^^ 

0<t'<[o— 1]+ t'\ 

To start the proof of (3) we write 



sup 

X 



I A(0'^) (t; xi, . . . , xf, X) I = r{t; xi, . . . , xj;X) = e^Ws{t) 
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and let a' = (t{xi, . . . , Xj] X) = (t{x, xi, . . . , xj; X), Then for r' e [0, [a' — 1]+] and 
t > we obtain, as in the proof of Proposition 2.7 

hence 

|Ar^)(t;xi,...,x,;X)| <M,e^-^N inf 

0<T<a' t! 

and the proof of (3) continues as the proof of (1). 

The case i > (taking now d{x, X) > i) is treated as in the proof of Proposition 
2.7. D 
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